home *** CD-ROM | disk | FTP | other *** search
/ Shareware Grab Bag / Shareware Grab Bag.iso / 004 / mathstat.arc / S16.BAS (.txt) < prev    next >
Encoding:
GW-BASIC  |  1983-11-14  |  896 b   |  42 lines
  1. 10  PRINT "GEOMETRIC REGRESSION"
  2. 20  PRINT
  3. 30  PRINT "NUMBER OF KNOWN POINTS";
  4. 40  INPUT N
  5. 50  J = 0
  6. 60  K = 0
  7. 70  L = 0
  8. 80  M = 0
  9. 90  R2 = 0
  10. 100  FOR I = 1  TO N
  11. 110  PRINT "X,Y OF POINT" ;I;
  12. 120  INPUT X,Y
  13. 130  Y = LOG (Y)
  14. 140  X = LOG (X)
  15. 150  J = J + X
  16. 160  K = K + Y
  17. 170  L = L + X ^ 2
  18. 180  M = M + Y ^ 2
  19. 190  R2 = R2 + X * Y
  20. 200  NEXT I
  21. 210  B = (N * R2 - K * J) / (N * L - J ^ 2)
  22. 220  A = (K - B * J) /N
  23. 230  PRINT
  24. 240  PRINT "F (X) ="; EXP (A); " * X ^";B
  25. 250  J = B * (R2 - J * K/N)
  26. 260  M = M - K ^ 2/N
  27. 270  K = M - J
  28. 280  PRINT
  29. 290  R2 = J/M
  30. 300  PRINT "COEFFICIENT OF DETERMINATION (R^2) ="; R2
  31. 310  PRINT "COEFFICIENT OF CORRELAION ="; SQR (R2)
  32. 320  PRINT "STANDARD ERROR OF ESTIMATE ="; SQR (K/ (N-2))
  33. 330  PRINT
  34. 340  PRINT "INTERPOLATION: (ENTER X=0 TO END PROGRAM)"
  35. 350  PRINT "X =";
  36. 360  INPUT X
  37. 370  IF X = 0  THEN 410
  38. 380  PRINT "Y ="; EXP (A) * X ^ B
  39. 390  PRINT
  40. 400  GOTO  350
  41. 410  RUN "STAT"
  42.